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Since its commissioning in autumn 2012, Tunka-Rex, the radio extension of the air-Cherenkov 
detector Tunka-133, performed three years of air shower measurements. Currently the detector 
consists of 44 antennas connected to air-Cherenkov and scintillator detectors, respectively, placed 
in the Tunka valley, Siberia. Triggered by these detectors, Tunka-Rex measures the radio signal up 
to EeV-scale air-showers. This configuration provides a unique possibility for cross-calibration 
between air-Cherenkov, radio and particle techniques. We present reconstruction methods for 
the energy and the shower maximum developed with CoREAS simulations, which allow for a 
precision competitive with the air-Cherenkov technique. We apply these methods to data acquired 
by Tunka-Rex in the first year which we use for cross-calibration, and we compare the results with 
the reconstruction of the energy and the shower maximum by Tunka-133, which provides also a 
reconstruction for the shower core used for the radio reconstruction. Our methods have shown that 
the atmospheric depth of the shower maximum (W max ) can be reconstructed with a precision of 
better than 40 g/cm 2 for high quality events, in some cases even when only three antenna stations 
have signal. The energy precision is comparable with the air-Cherenkov precision of 15%. Soon 
the results will be checked with the independent data of the second year. 
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1. Introduction 

Radio detection of high-energy cosmic rays recently has become a competitive technique. 
Provided a simple external trigger it reaches a duty cycle around-the-clock. A first generation 
of modem radio experiments [1,2] has opened the prospects cosmic-ray detection with energies 
higher than 0.1 EeV. The dense radio-astronomical observatory LOFAR demonstrated good agree¬ 
ment between the amplitude footprint predicted by CoREAS and real measurements [3]. After the 
confirmation that the radio signal is sensitive to the air-shower parameters energy and shower max¬ 
imum, the main remaining question is the applicability to large scales the radio technique for the 
detection of ultra-high energy cosmic rays. Currently there are two experiments aiming to answer 
this question: AERA [4] and Tunka-Rex [5]. 

Tunka-Rex is the radio extension of the air-Cherenkov detector Tunka-133 [6] placed in Siberia, 
near Lake Baikal. The host experiment Tunka-133 is successfully taking data since 2009. It op¬ 
erates during winter in moonless nights, reaching a total duty cycle of less than 500 h/year. The 
detector array consists of 19 clusters placed on 1 km 2 area, each equipped with 7 photomultiplier 
tubes (PMT) with spacing 80 m. Later it was extended by 6 externals clusters increasing the effec¬ 
tive detection area up to 3 km 2 . Recently an additional array of scintillators, Tunka-Grande, and 
air-Cherenkov extensions under the name of TAIGA have been installed [7]. 

The radio extension (Tunka-Rex) consists of 44 antennas, 25 of them are connected to Tunka- 
133 cluster centers. Another 19 are connected to the scintillator extension currently under com¬ 
missioning. The Tunka-Rex detector works in slave mode, i.e. it receives trigger from the air- 
Cherenkov array; all events are recorded in parallel. This provides the ultimate possibility of 
cross-calibrating the radio and the air-Cherenkov techniques. Since Tunka-133 has a high resolu¬ 
tion (15% for the energy and 28 g/cm 2 for the shower maximum reconstruction), Tunka-Rex has 
the potential to determine the precision of the radio technique. 

Tunka-Rex started its data acquisition in October 2012. Until now, three seasons of data 
acquisition were completed (2012/2013, 2013/2014 and 2014/2015). Our collaboration decided to 
follow a semi-blind analysis paradigm: energy and shower maximum reconstruction by Tunka-133 
is blinded starting from the second season. In the present work we show the results of the first 
season used for tuning and testing the reconstruction method. These methods are used to predict 
the energy and shower maximum for the second season enabling an independent cross-check. 

2. Reconstruction of air-shower parameters 

For the reconstruction of air-shower parameters we use a common method exploiting the lat¬ 
eral distribution of amplitudes. In principle, if one knows the shower axis and core position, this 
method requires only the amplitude at two antenna stations for the energy and shower maximum 
reconstruction, since mainly these two shower parameters are defining the lateral distribution. The 
amplitude is sensitive to the energy, and the shape of the lateral distribution to the shower maxi¬ 
mum (see, for example, Ref. [8]). For the radio detection, the main difficulty rises from two features 
of the generation of radio emission. First, the interference of the geomagnetic effect [9] and the 
charge excess [10] introduces an azimuthal asymmetry in the lateral distribution: the distribution 
becomes two-dimensional. The second effect comes from the features of the signal propagation in 
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a medium (atmosphere) with refractive index n r > 1 : the lateral distribution contains Cherenkov- 
like features. Taking these features into account, we developed an appropriate methods to describe 
measured lateral distributions. Our study is based on CoREAS [11] simulations and described in 
Ref. [12]. 

2.1 Asymmetry 

In Ref. [3], it was shown that the charge excess contribution to the radio emission leads to 
significant complication of lateral distribution. Nevertheless, this distribution can be fitted with 
a two-dimensional function [13]. This solution is very general, but hardly applicable to sparse 
detectors operating near threshold, i.e. a typical event contains only 3-5 antenna stations with 
signal. We provide a simpler approximation of the lateral distribution and parameterize it with a 
one-dimensional function reducing the number of free parameters. 

The squared signal amplitude at a certain distance has the form 

= S’q ((sin a g + £ cos </> g ) 2 + £ 2 sin 2 0 g ), (2.1) 

where £ is the asymmetry defined as fraction of strengths of the Askaryan relative to the geomag¬ 
netic, a g is the geomagnetic angle (angle between vector of magnetic field and shower axis), is 
the azimuth of an antenna station in the shower plane. 

As we can see, the one-dimensional LDF transforms itself to a two-dimensional one when tak¬ 
ing into account the contribution from the charge excess phenomena: £ > 0. To reduce the number 
of dimensions back to one, we define a special operator K eliminating the azimuthal dependence 

K = \Je 2 + 2£ cos (j > g sin a g + sin 2 a g , (2.2) 

The remaining question is the value of the asymmetry £, which should be used for the correction. 
In Ref. [12] it was shown with CoREAS simulations that a constant value £ = 0.085 is sufficient 
for the correction of the lateral distribution function for the situation of Tunka-Rex. 

2.2 Lateral distribution function (LDF) 

As mentioned, due to the refractive index a Cherenkov ring is present in the lateral distribu¬ 
tion [16, 17]. Thus, we use the simplest function which can describe the Cherenkov ring [18] 

Sir) = £’r 0 exp(ai(r-ro)+a 2 (r-ro) 2 ), (2.3) 

where r is the distance to shower axis. A problem rising from this parameterization is connected 
to the strong correlation between the parameters a\ and a 2 . When the number of fitted points is 
close to the number of fit parameters, in addition to large error bars, the fit can converge to false 
minima, or give large uncertainties for parameters with strong correlation. Unfortunately, most of 
the Tunka-Rex events satisfy these conditions. For this reason we decided to reduce the number of 
free parameters in the LDF by fixing the parameter a 2 to a value depending on zenith and energy 

a 2 (8,E pT ) = fl 2 ()(Fpr) + 021 (Fpr) cos 0 > (2.4) 

«20 = «200 + 0201 Epr ■, d 2 \ = ^210 + 0211-Epr j (2-5) 
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with «200 = 0.19 • 1(U 5 m~ 2 , ^201 = —1.63 ■ 10 -5 m~ 2 /EeV, 0210 = —3.45 • 1(U 5 m~ 2 , ai\\ = 
2.44 ■ 10 5 m 2 /EcV obtained by the CoREAS simulations. The primary energy E pr for this pa¬ 
rameterization is estimated with a simple exponential LDF after correction for the asymmetry. The 
final energy estimator is described in the next section. 

2.3 Energy and shower maximum reconstruction 

The detailed explanation of formulas used for the reconstruction of air-shower parameters is 
given in Ref. [12]: here we only give the parameters obtained from the model. The energy can be 
reconstructed by probing the signal amplitude at a defined distance r e 

E VI = K L g{r e ), (2.6) 

where parameters have the following values: r e = 120 m, = 884 EeV/(V/m) 1 . 

For X max reconstruction we use the parameterization 

X max =X det /cos6-(A + B\og(ai(r x ) + b)), (2.7) 

where A c i el = 955 g/cm 2 is the atmospheric depth of the detector. The parameters have the following 
values: r x = 180 m, A = —1864 g/cm 2 , B = —566 g/cm 2 , b = 0.005 m -1 , 

3. Event reconstruction and comparison with Tunka-133 

We perform the cross-calibration for the events of 2012/2013, and make the prediction for 
the data of the season 2013/2014, which is still blinded. After unblinding we will compare the 
prediction given by Tunka-Rex and Tunka-133, and present the conclusion for the precision of the 
air-shower reconstruction. These results will be published in another paper soon. In the present 
work we present results of the first season in comparison with Tunka-133 reconstruction for the 
energy and shower maximum. 

3.1 Data acquisition and signal selection 

All events are triggered by air-Cherenkov detector, the radio detector is read out in parallel. 
As main software for data analysis we use a modified version of the radio extension of the Auger 
Offline framework [14], For the reconstruction of the radio signal we apply antenna patterns, and 
hardware responses based on calibration measurements [15], As result, we obtain signal traces for 
each antenna station in the units of the electrical field. We require at least three stations with signals 
passing a signal-to-noise (SNR) ratio cut: SNR = S 2 /N 2 > 10, where 5 is a signal amplitude and N 
is the RMS of noise calculated in a noise window. With simulations we parameterized the impact 
of the noise on the measured amplitude and correct for this. 

Using the time information of these stations we can reconstruct the air-shower arrival direction 
and compare it with the reconstruction given by Tunka-133. Events with big deviation (> 5°) in 
arrival reconstruction are considered as false and rejected. As last step we remove outliers from 
the lateral distribution and fit this distribution with our parameterization, if it still has at least three 
antenna stations with signal. An example event is shown in Fig. 1. 

'Here we use Ky which differs from the value K = 422 EeV/(V/m) in Ref. [12]. This is because we simplified the 
equation of Ref. [12] by setting the exponent b = 1 instead of 0.93. 
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Figure 1: Example event reconstructed by Tunka-Rex. Left: Footprint of the air-shower. The color code 
shows the arrival time. Points show clusters without antenna in 2012/2013. Right: Lateral distribution with 
fit. Grey points show the antennas not passing the SNR cut, the dotted point is a rejected antenna with false 
signal. 

In the season 2012/2013, the Tunka facility had effectively 280 h of measurement time. The 
full reconstruction of Tunka-133 is available only for events with zenith angles 6 < 50° due to 
design features. Thus, for the Tunka-Rex cross-calibration we used only these events. After the 
cuts described above we have 91 events. 

3.2 Energy reconstruction 

We present the results of the energy reconstruction given by formula (2.6). The spread between 
values given by Tunka-Rex and Tunka-133 is 19 ± 3% when using all 91 events (see Fig. 2). 

3.3 Shower maximum reconstruction 

The reconstruction of the shower maximum is connected to the slope of the lateral distribution. 
Dense detectors have shown good sensitivity to the shower maximum [3, 13], but the reconstruction 
with sparse detectors as Tunka-Rex is more complicated. Events with small number of antenna 
stations can be used to reconstruct the shower maximum, only if the fitted slope (parameter a\ of 
Eq. (2.3)) has a small uncertainty. However, a small number of antenna stations means that the 
lateral distribution has points close to each other, which leads to an underestimation of the slope, 
i.e. an overestimation of the distance to the shower maximum (see dotted events in Fig. 3). Since 
we probe the slope at the distance of 180 m, we exclude all events which have not antennas at 
a distance > 200 m to the shower axis, thus cut removing 27 events. The fitting uncertainty of 
a i propagates to the uncertainty of the shower maximum reconstruction <7x max . We define a high 
quality cut for this uncertainty <7x max < 0$^ . For the present analysis we set = 50 g/cm 2 . 
After this cut 25 events survive (bold points in Fig. 3). The mean deviation between Tunka-133 and 
Tunka-Rex for this quality data set is 51 ±7 g/cm 2 . Decreasing the cut threshold crj^ we further 
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Radio energy (EeV) Deviation 

Figure 2: Correlation between the energy reconstructed by Tunka-Rex and by Tunka-133, and histogram of 
relative deviation. The deviation is calculated as energy difference divided by average. 

increase the precision of Tunka-Rex losing more events. The best result reached with our statistics 
is 36 ±7 g/cm 2 for 15 events. This means that the uncertainty is consistently estimated. If 
the Tunka-Rex uncertainty for X max would be equal to the Tunka-133 uncertainty we would expect 
a mean deviation of y/2 ■ 28 « 40 g/cm 2 . 

4. Conclusion 

Tunka-Rex is an excellent experiment to provide a benchmark for the radio detection of cosmic 
rays. The detector is a sparse array of antennas, and operates near the detection threshold. Most 
events have low SNR and contain only about 3-5 antennas with signal. On the other side, the 
experiment is hosted by a precise air-Cherenkov array, which reduces the systematic uncertainty 
for cross-calibration. 

Our study has shown that taking into account all significant effects of the radio emission allows 
to work under these strong restrictions. One of the most important questions in this topic is the 
precision for the reconstruction of the shower maximum. We have shown that a sparse array with 
spacing of 200 m is able to reconstruct the shower maximum with uncertainty better than 40 g/cm 2 
even for some events containing only three antennas, if the event contains antennas with signal 
further than 200 m from the shower axis for a high-quality selection. Tunka-Rex within its energy 
range and detector size has uncertainty of about 50 g/cm 2 after cross-calibration with Tunka-133, 
which has a precision of 28 g/cm 2 . Taking all the aforesaid into consideration, one can state, that 
radio detection is economically feasible for ultra-high energy cosmic rays detection. 

Tunka-Rex recently has been extended with additional 19 antennas connected to Tunka-Grande, 
the scintillator extension of the Tunka facility. Triggering by the particle detectors will increase the 
duty cycle by an order of magnitude. This configuration provides a possibility for cross-calibration 
between three air-shower detection techniques, as well as an increase of the statistics of the Tunka 
facility in the energy range around 1 EeV. 
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Figure 3: Correlation between the distance to the shower maximum reconstructed by Tunka-Rex and by 
Tunka-133m and histogram of the relative deviation. Both reconstructions are in agreement within the 
uncertainties. The mean deviation is 51 ± 7 for the high and 130 ± 12 g/cm 2 for the low quality cuts. The 
deviation is calculated as A max difference. 
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